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Synopsis 

We solve massive gravity field equations for all locally homogenous and vanishing scalar invariant 
(VSI) Lorentzian spacetimes, which in three dimensions up to fibering and warping are the building 
blocks of constant scalar invariant (CSI) spacetimes. At first, we provide an exhaustive list of all 
Lorentzian 3D homogeneous spaces and then we determine the Petrov type of the relevant curvature 
tensors. Among these geometries we determine for which values of their structure constants they are 
solutions of the field equations of massive gravity theories with cosmological constant. The homogeneous 
solutions founded arc of all various Petrov types : Ic, Iw, II, HI, Dt, Ds, ^, O; the VSI geometries 
which we found are of Petrov type ///. 
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1 Introduction 

ft is known that, in three dimensions, Einstein gravity theory does not possess any local physical degrees 
of freedom. However Deser, Jackiw and Templeton proposed a modification of the theory consisting 
to add to the usual Einstein-Hilbert Lagrangian a Chern-Simons term |2l [3] , defining so what is known 
as topological massive gravity (TMG). So they provided a chiral three-dimensional gravity theory that 
acquires massive spin-2 excitations and mediates finite-range interactions. This modification also proved 
to be an ideal example for a consistent quantum theory of gravity: it is unitary (4] [5] and super- 
renormalisable 0. 
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On the other hand BergshoefF, Hohm and Townsend [71 [H] have recently obtained another type of 
massive gravity theory (NMG): a parity preserving theory that describes (on a Minkowski background) 
the propagation of a massive positive energy spin-2 field, but now of both helicities ±2. This theory 
possesses all the virtues of the TMG, so it may be considered as another consistent theory of quantum 
gravity. They have also considered the "merging" of both theories, producing a general massive gravity 
theory (GMG) which involve two spin-2 helicity states with different masses (parity-violating) and, as 
the previous ones, can be also extended by adding a cosmological constant. 

Our work is motivated by the paper of Chow et al. 9J which provided a review of a large set 
of solutions of topological massive gravity (with cosmological constant). In their paper these authors 
provided a 3D variant of Petrov classification and showed that all the solutions they founded in the 
literature at this time are of Petrov type D or N , corresponding to locally squashed AdS^ or AdS pp- 
wave metrics. Moreover they also proved that all Petrov type D solution of TMG actually are biaxially 
squashed AdS^ metrics. Moreover, in a companion paper jlOj these authors obtained new solutions of 
the topologically massive gravity equations by considering Kundt metrics jTH [T2] [see also Chapters 28 
and 31 of [13 ]. These metric are characterised by the existence of a null geodesic vector field (that in 
3D is automatically shear- free and twist-free). The TMG solutions belonging to this class of metric are 
generically of Petrov type //, but also in some special cases of Petrov types I?, ///, N and O. The 
classification of the homogeneous solutions to the classical vacuum equations of TMG were found by 
Ortiz in |14j and was recently generalised for non- vanishing cosmological constant |15j . 

On the other hand, it was proved in |16) that a Lorentzian 31? spacetime which is locally a constant 
curvature space (CSI) can be constructed, by means of fibering and warping, from locally homogeneous 
spaces and a subclass of CSI which are the vanishing scalar invariant spaces (VSI). This motivated our 
restriction to the geometries we considered in this work. 

The paper is organised as follows; In section [2] we summarise the various field equations of the 
massive gravity theories. We start section [3] by revisiting all the three dimensional homogeneous spaces 
with Lorentzian signature and determine the Petrov type of the relevant curvature tensors (listed in the 
appendix), which will prove to be useful for solving the equations and identifying the solutions. We 
then proceed in section |4] to solve the equations of massive gravity theories for homogeneous spaces with 
cosmological constant. In the last section [5] we consider solutions of massive gravity theories for VSI 
geometries. 
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2 Massive Gravity Theories 



Let us first remind the Einstein-Hilbert action, 

Seh = ^ J V\9\{R~'^^)'i^x , (2.1) 

where A is the cosmological constant and k is the gravitational couphng with mass dimension [k] ^ — \ ■ 
In what follows we adopt the mostly plus expression of the metric, define the curvature tensors so that the 
curvature of the Euclidean round sphere equipped with its positive definite metric as positive curvatur^ 
and fix the spacetime orientation by choosing eoi2 = +1- 

2.1 Topologically massive gravity 

Topological massive gravity is obtained by adding a gravitational Chern-Simons term to the Einstein- 
Hilbert action 

Stmg := Seh H ^ Scs , (2-2) 

which is expressed through Christoffel symbols of the spacetime metric g while is a new coupling 

constant with mass dimension one. 

The classical equations of motion read as 

Rpu - ^gt^uR + ^gt.u + ^ C^^ = , (2.3) 



where is the Cotton- York tensor : a symmetric, traceless and divergenceless tensor. An immediate 
consequence of the equations of motion is that the traceless part of the Ricci tensor and the Cotton- York 
are proportional 

Si_iv H Cfj_y — , Sfj^i, :— Rfiu — ^g^yR , (2-4) 
^ o 

and accordingly of the same Petrov type. 
2.2 New massive gravity 

This theory is defined by adding a quadratic curvature term [7] to the Einstein-Hilbert action 



Snmg '■= Seh — -. — ^ Sqc , (2.5) 
4 k 



Sqc = J V\g\ {r^.uR^" - ^i?') d^x 



^ In others words according to the conventions : R'^^ p a = dpV'^ + ■ ■ ■ \ Ru a = R''v p c 



where ^ is a coupling constant with mass dimension two. The corresponding equations of motion read 
Rt.u - ^.9m'v-R + A.g^^ -^K^^^O , (2.6) 

where K^^, is a symmetric and divergenceless tensor whose trace coinsides with Lagrangian density of 



Sqc ( [23| ), namely 

K = g^'-K,, = i?^,/?^"- - . (2.7) 

A consequence of the equations of motion is that traceless parts of the Ricci tensor and i^^^ tensor 
are proportional 

Sf_tu = 2^^t^^ ' -^A'l' ^t^f ^ 7^9fivK . (2.8) 

Thus, the Petrov type of the traceless part of the Ricci tensor and of K^^, have to coincide in the 
framework of the new massive gravity. 

2.3 General Massive Gravity 

This is defined by the merge of topological and new massive gravity theories 



Sgmg Seh H Scs — Sqc , (2.9) 



and the corresponding field equations read 



R^.u-^9^.uR + -f^9^,u + ^Cf,^- ^Kf,^ = . (2.10) 

To understand the field content of the latter equations, we have to linearise them around a fixed back- 
ground. This was performed around Minkowsky and AdS:^ backgrounds in [7] and |17| respectively. In 
particular, it was shown in the latter that there are two massive spin-2 modes, which are stable if 

e(^ + V)^2AAi2 . (2.11) 

Thus, for pure NMG, when /j, goes to infinity, AdS^ is stable if ^ ^ ^. In what follows we shall not 
restrict the sign of the coupling constants, but just discuss them according to the parameters appearing 
in the expressions of the metrics we obtain. 

3 Homogeneous geometries 

The strategy we adopt to obtain all homogeneous spaces, solutions of the massive gravity field equations 
was the one advocated a long time ago by Ozsvath [TH]. We take as metric components the canonical 



form 

(r?„/3)=diag.(-l,l,l) (3.1) 
but start from arbitrary structure constants of the Lie algebra and fix them using only /so(l, 2) transfo- 



rmations that preserve the expression (3.1 ) of the metric. Our strategy, is to first review the classification 



so to use it in the field equations and then solve the resulting algebraic equations. 
3.1 Bianchi classification revisited 

As it is well known (see for example [T^l) the Lie algebra of the (right) invariant fields 

[Uip]^-Cl^^, , (3.2) 

whose structure constants satisfy the Jacobi identity Cf-^C'^^^e"'^'^ = 0, allocates in two classes. The 
unimodular Lie algebras, such that C^^ — and the non-unimodular ones such that \C'^^ — ka ^ 0- 
For the three dimensional real algebras (first classified by Bianchi), their structure constants can be 
parametrised in terms of the vector components ka and a symmetric tensor densitjj^ n" ^ (see |13) and 
refs therein). 

C^^=sp^cn'^" + k06^-k^S'§ (3.3) 

while the Jacobi identity reduces to the condition k^n"'^ = 0. Thus, the /so(l,2) classification of the 
structure constants reduces to the classification of symmetric tensor densities, a problem solved in [19] 
(see also refs [20l[2T]), annihilating a vector. 

Chow et al. [9! suggested to classify solutions of TMG according to the Segre classification of the 
traceless Ricci tensor Sa p , which is equivalent to the Petrov classification of the Cotton- York tensor 
but not necessarily in the frame work of NMG. Here after we present the various canonical forms of the 
Lie algebra obtained using /so(l,2) transformations and the Segre-Petrov type of their traceless Ricci 
and Cotton- York tensors. On a practical level, to make this classification it is not necessarily to know 
exactly the eigenvalues of the tensor we want to classify when they are all different. So we have just 
compute the discriminant of the characteristic equation. If it is positive, we know that two eigenvalues 
are complex conjugate and one is real and so corresponds to the case Ic (see ref.[9] for the notations); if 
it is negative the three eigenvalues are real and distinct, corresponding to the case Jr. It is only when it 
is zero, in which case we know that at least one eigenvalue is double, that further analysis is required to 
determine the Jordan form of the tensor. But the degeneracy of the eigenvalues greatly facilitates the 
analysis. In this case if the characteristic polynomial reduces to its cubic term P{X) = A'^, the tensor 
is of Petrov type O, N , or ///; otherwise when the roots are equal, but non vanishing, the tensor is of 
Petrov type Dg, Dt or //. 



In what follows we shall call them respectively structure vector and structure tensor density. 
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3.2 Unimodular Lie algebra 



Four types of normal forms are possible 
• Type / 



K^) = I 6 c I , (3.4) 
\0 c b J 

this form is equivalent to diag. (a, 6 + c, & — c) , but the latter appears to be more easy to handle for 

solving the field equations. By changing the orientation (in which case fi also changes its sign), we 

always may assume that a > if it is non zero, otherwise that & > if non zero, etc. 

Obviously this structure tensor density may correspond to any of the unimodular Bianchi types 

(/, //, VIq, VIIq, VIII, IX) of real Lie algebras (see for instance Ortiz [Hj). 

The matrix component of Ricci tensor in the invariant frame, of the metric so defined reads 

-2-402 

(^a0) = ^| a(a + 2&) -2(a + 2&)c I , (3.5) 





-2(a + 2&)c a{a + 2b) 
and its trace is equal to i? = ^ (a^ + iab + Ac^). Accordingly the traceless part of this tensor is 



(Sa/i) 



( I (a2 + a6-2c2) 

i(a2 + a6-2c2) -(a + 26)c I . (3.6) 



\ -(a + 26)c 

To determine the Petrov type of this tensor, we have to obtain the Jordan form of the matrix 
of components S^. Fortunately, we have not to handle explicit solutions of the third degree 
characteristic polynomial 

P{X):^dei[\61-S'^] (3.7) 

in the general case. The traceless condition implies that this cubic polynomial, always will be of 
the form : P{\) = \^ + p\ + q. Its discriminant, which is defined as A := ^ + 1^, partially fixes 
the number and the nature of the different roots of P{X) — {). If A > one root is real and the 
two others complex conjugate; if A < the three roots are real and distinct; if A = 0, at least 
two roots are equal. Accordingly, when A > the traceless tensor will be of Petrov type Ic and 
when A < of Petrov type Jr. When A 0, the tensor is of special Petrov type. It is of type // 
or D when p 7^ 0, and of type ///, or O when P{\) — . Its precise determination will need 
more investigation, but things are greatly facilitated as at least one of the eigenvalue of the tensor 



is degenerate. For example, in case of the structure tensor density (3.6 1, we obtain 



As = (a + 2 6)2 c2 ((a + bf - c^f {a' - Ac^ , (3.8) 
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which generically is negative and thus define a tensor of type /«; exceptions occur when it vanishes. 

The same discussion can to be preformed on the Cotton- York tensor 

-a2(a + 6)-46c2 \ 



\ (a3 + 6a2+46c2) c(2c2 - ia^ + 46^ + 2a6) 
c(2c2- ia2+462 + 2a6) \ {a^ + ba"^ + Abc^) ) 



(3.9) 



The results of this analysis, both for the Sa p with Ca p and Ka p tensors, are summarised in tables 
1-6. Thus, solutions of TMG or NMG can be easily found by comparing the Petrov classifications 
of Safi with or Kafi respectively. However, this is not a priori the case for GMG. 



Type II 

V + a V 
V I/ — aOl , u = ±1 
6 

which does not admit timelike eigenvector but a double null vector. 



K/) = 



(3.10) 



Diagonalising n/j with a GL{3, M) transformation we see that if a = 6 = 0, this corresponds to a 
Bianchi // type of Lie algebra, if 6 = but a ^ or a = and ub < to Bianchi VIq, if a = 
and > to Bianchi VIIq, and otherwise to Bianchi VIII. 



The Ricci and Cotton- York tensors road 

-{2a + b){b-2v) -2{2a + b)v 

-2{2a + b)v {2a + b){b + 2u) 

-9 



{Rap) = 2 



(3.11) 



Type /// 



{{8a'^ + 4:ba~b'^)u + b^{a + b)) {Sa'^ + 4ba - h^) u 

-{8a'^ +4ba-b'^)iy {8v - b {b - 4i^) a - b^b + i^)) 

26^(0 + 6) 



/ a 1 \ 

1 -a 1 
^ 1 -a J 



(3.12) 



(3.13) 



which admits a triple null vector. 

This structure constant density may only corresponds to Lie algebra of Bianchi types VIo if a = 
and for o 5^ it is VIII. 



The Ricci and Cotton- York tensors are 



( 



\ 



2-2 


—a 


2 


-a 




a 


2 


a 




6 a 




—6 a 


2"' 







6 a 




—6 a 



\ 



(3.14) 



(3.15) 
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Type IV 



/ v 0\ 



(n'^P) = „ a , (3.16) 

^ 6 y 

where < 4i/^ and which only has one simple spacelike eigenvector, but no timelike or null 
eigenvector. 

Here again, only the structure constants of Lie algebras of Bianchi types VIq if 6 = and for 6 ^ 
it is VIII. The corresponding Ricci and Cotton- York tensors read 



-{a-bf 2{a-b)v 
{R^p) = ^\ 2{a-b)v b'^-a^ 



1 



(3.17) 



(C7«) = I -„ (2(a2 _ ,,2) _ 1 (a + bf) a{a^ - 2i.'') - i6(a' + &') 



(3.18) 



3.3 Non-unimodulcir Lie algebra 

Here we have to consider three possibilities: the vector of components ka is timelike, spacelike or null; 
four type of normal forms will occur. 



Timelike ka ■ We choose the frame such that k^ = {k, 0, 0). The Jacobi identity implies that ^ 
is a spacelike symmetric tensor that can be diagonalised by rotation in the [1, 2] plane. But this 
normal form is not the most suitable for the resolution of the field equations and we prefer to use 
the following form 



/ \ 



(3.19) 



a b 
\ b a J 

Obviously, all non-unimodular Lie algebra may lead to this form of the stnicturc tensor density. 
More precisely we have a Lie algebra of Bianchi type V if a = 6 = 0, of type IV if a = ±6 ^ 0, 
and otherwise of types ///, Vlh for |6| > \a\ or Vllh for \a\ > \b\ with h = k'^ /{a^ — 6^). 



The corresponding Ricci and Cotton- York tensors read 

-(62 + P) \ 

{Rap) = 2\ k{k + b) -ab 

-ab k{k-b) j 



( -2ab'' 



iCf) = 2 











a6(&-3fc) b(2c? + -k'^) 

\ 6(2o2 + 62- a6(6 + 3fc) / 



(3.20) 



(3.21) 
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Spacelike ka : We choose the frame such that ka = (0, 0, k). Then the structure tensor density may 
take three different canonical forms according to that it admits a timelike (and thus a spacelike) 
eigenvector 



/ a \ 



(3.22) 



6 
\ / 

corresponding to any type B Bianchi space, namely: Bianchi type V if a = b = 0, type IV if a or 
b non-zero, and otherwise type ///, Vlh for o6 < or Vllh for ab > with h = fc^/ (a b). If it has 
a double null eigenvector : 



/ ly + a 



inrn) = 



\ 



(3.23) 



u u — a 
\ / 

corresponding to Lie algebra of Bianchi type IV if o = 0, otherwise to Bianchi type /// and Vlh 
with h = —k^/a^. If the structure tensor density did not have any other eigenvector, it can be put 
in the form : 



/ \ 



infill) = 



(3.24) 



1/ o 

\ y 

and corresponds to Bianchi types /// and Vlh with h = —k'^/u^. In SI case, the Ricci and 
Cotton-York tensors read 



1 



a" - 6" + 4fc^ 



'2k (a + b) 



(-Ra/3) = ^| -2k{a + b) a'^-b'^-4k'^ 



\ 



(3.25) 



{a + bf-4P J 



^/3y 2 



in case 5/7 



-(a + &)(2a^-a& + 62 + 2fc2) ?,k(c?-\?) 

-3^(0^-6^) (a + 6)(a'-a6 + 2(6" + fc^)) 

(o + 6)(o2-62) 



k^ + (k-\'a)v -iy{a + k) 
(Raff) =2 \ -v{a + k) u{a + k)-k'^ 




(3.26) 



(3.27) 



and in case SIII 



1 -1 

-2v{a + k){2a + k) | 1 -1 




(3.28) 



{R. 



■afi) 




a{v — k) 



1„2 



2k{k-u) 







2(fc2 + Z/2) 



(3.29) 
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( -a{\a? + k{k-Zv)) - z/) (4i/(A; + - Sa^) 

-\{k-v){Av{k + v)-'ia'^) a{a^ + P -Ski^-2iy'^) 

a{2u'^-\a 



(3.30) 



Lightlike ka ■ Here, without lost of generality, we may assume ka = (1, 1,0). Using the Jacobi 
identity we obtain the expression 



a 


—a 


b 


—a 


a 


-b 


b 


-b 


c 



(3.31) 



corresponding to Bianchi type Fifo = 6 = c = 0, type IV if 6^ = ac and Bianchi type ///, Vlh for 
b"^ > ac or Vllh for 6^ < ac with h = l/{ac — b^) in the other cases. Without lost of generality, it can 
still be simplified, by performing an appropriate null rotation around ka, that leads to 



fa -a 0\ 

—a a or 
c y 

The Ricci and Cotton- York tensors are 

/ c(a - ic) - 26(1 + b) 



/ b 

0-6 
^ b -b 



{Rap) 



\ 



ac-2b{l + b) {l + b)c 
ac- 26(1 + 6) c(a+ ic) -26(1 + 6) (1 + 6)c 
(1 + 6)c (1 + 6)c 



2^ 



(3.32) 



(3.33) 



-c (262 + 36+ 1 + ic(c- 
c(262 + 36+ 1 - iac) 



a)) -c(262 + 36+l- iac) f(l + 6)c2 



c (262 +36 + 1 - ic(c + a)) -|(l + 6)c2 



(3.34) 



-|(l + 6)c2 -i(l + 6)c2 

4 Solutions of the field equations 
4.1 Simply transitive groups 

In the framework of the homogeneous spaces, whose group of isometrics is strictly transitive (Bianchi 
spaces), the field equations reduce to algebraic equations. Our strategy in solving the equations, is 
to first obtain the link between the cosmological and structure constants A and fi {resp. ^ ) and the 
structure constant parameters and then inserting them into the field equations and discuss the remaining 
constraints that have to be satisfied. 

We shall provide some details in the first case, whereas in the other cases we shall just present the 
solutions of the equations. 



Unimodular Lie algebra 
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Type / 

★ TMG: Using the expression (3.4 1 of the tensor density defining the structure constants, 
the TMG field equations reduces to three independent equations 



2 (a + 6) + 8 6 + (a^ - 4 + 4A)^ = 

— +ab + c^-3A = 

c(-a2 + 4a 6 + 8 6^ + 4 - 2(a + 26)^) = 



(4.1) 
(4.2) 
(4.3) 



which are easy to solve. 



First let us assume c = 0, so that eq. (4.3) is trivially satisfied. We deduce from the other two 
that: 

1 



A= ^a(a + 46) 



from which we obtain the solution of Petrov type Dt 



^2 -27 A 
6 /i 



c = 



(4.4) 
(4.5) 

(4.6) 



To isolate the value (4.51 of ii in eq. (4.1) we have assumed that a + 4A 7^ 0. If a + 4A 0, 
the parameter /i remains undetermined and the field equations are satisfied if the cosmological 
constant is negative and 



-b = 2V-A 



(4.7) 



or if A — 0, for a = and b arbitrary. Of course these last two solutions are conformally flat 
(Petrov type O). More precisely if a = the solution is flat; otherwise it is AdSs ii a = —b ^ 0. 



If we assume c 7^ 0, the values of A and ^ provided by the first two field equations (4.1), (4. 2) 
become 

1 



^-12'" 



2 +4a6 + 4c2) 

0^(0 + 6) + 46c2 
a{a + 6) - 2c2 



Inserting these values into the third non trivial field equation (4.3) we obtain 

((a + 5)2-c2)(a2+4afe + 4c2) =0 



(4.8) 
(4.9) 

(4.10) 



as we assume of course that /i 7^ 0. This leads to two solutions (changing the sign of c correspond 
to a reflexion of 6*^ or 9^ and thus to change the sign of /i) of Petrov type Dg'- 



27A-/i2 
6 fi 



b = 



5^*2 - 27 A 
12ji 



±{a + b)=±'- 



9 A 



(4.11) 
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and a third one, of Petrov type /r, non flat, but with vanishing cosmological constant and 

b=^{fi + a) , = -^a{3a + 2^i) , (4.12) 

which restrict the parameter a to belongs to the intervaQ a e [0, — • 

* NMG: The strategy is the same, but the equations a little bit more cumbersome. The 
field equations lead to 

-210"* - 200^6 + SOafec^ + 256&2c2 + 80c^ + 8$ (a^ + 2a6 - 4A) = , (4.13) 
-eSa"* - SOaH + Sa^ (-2 6^ + 3 + 2 ^ +16 (l6 &^ + 5 - 4 ^ +4Af ) =0(4.14) 
c (5a3-2a2 6 + 40a62-8C (a + 2&) + 20ac2 + 6463 + 1046c2) =0 . (4.15) 

Here also three different types of solutions occur. 
If we assume c — then we obtain 

a{21a^ + 72 ab+ 16 b^) 
^ = 8(21a + 45) ' ^'-''^ 

^=^(21a + 4&) , (4.17) 

which leads to a solution of Petrov type Df : 

16 



"^-21 



(SC ±V3Ve (7A + 50) , (4.18) 

_ 42 AC ±73(21 A- 170 (7A + 50-66C^ 
^ - 4(A-e) ■ ^ ' 

We find more convenient to discuss this solution by introducing the parameter x = 6/a in terms 



of which we may rewrite eqs (4.161, (4.171 as 



21 + 72 a: + 16a;2 a . ^ . . 

A= , C=— (21+4a;) . 4.20) 

8 21 + 4x ' ^ 8 ^ ^ ^ ^ 

In order to have ^ > we need, in addition to a 7^ 0, that x > —21/4, which implies that we 
always have A > if -(2^/T5 + 9)/4 > x > -21/4 or a; > (2^15 - 9)/4 . 

If we did not assume c = 0, the first two equations ( 4.13^ , ( |4.14[ ) lead to much more complicated 
expressions 

^ _ 2ac2 (I5a3 + 28a2&- 128 53) +a^{a + b) {21a^ + 72 ab + 16b^) - 166c4(15a + 326) - leOc^ 
~ 8 (a2(a + &)(21a + 46) - 2c^ {3a^ + lOab + 6462) _ 40c4) 

(4.21) 

a^(a + 6) (21a + 46) - 2c2 (Sa^ + 10a6 + 646^) - 40c* 
^ 16 (a(a + 6)- 2 c2) ' ^ ' 

^When we write an interval as [a;o,a;i] we did not assume xg xi, but consider the union of the sets {a'|xo ^ a; 
xi} U ^ 2; ^5 xq} even, unless xo = xi, one of these two sets is always empty. This convention allows to avoid tedious 

(but elementary) discussion about signs. 
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inserting these in (4.151 field equation we find 

a{a + bf{a^ + 2ab- 46^) - (a^ + IQa^b + 12 ab^ + 8 b^)c^8 be* = (4.23) 

wfiose solutions arc = {a+b)^ = a^ll+x)^ and = a{a'^+2ab-4:b'^)/8b = a^(l+2a:-4x^)/8a;. 
From the first solution, which is Petrov type D^, we obtain 

which satisfy the requirement ^ > if a; ^ [—25/42,-1/2]. We can see immediately that this 
restriction is compatible with both signs of the cosmological constant : A < if x e] — (8\/T5 + 
57)/42, (8Vl5 - 57)/42[U] - 25/42, -l/2[, A = at the boundary of these intervals, excepted at 
X = -25/42 where it diverges; otherwise A > if a; ^ [-(8Vl5 + 57) /42, (8715 - 57)/42] U] - 
25/42,-1/2]. Thus here A is bounded from below (by approximatively — 0.3a^), but it can be 
chosen arbitrarily positive. 

For the second solution — a^(l + 2x — 4a;^)/8a; things are a little bit more involved. We obtain 

^ ^-^{Mx^~Ux' + 36x + 5) , A=|^— ^4^?47^ ■ (4-25) 

32a;^ ' 16x (64a;3 - 44x2 + 36a:; + 5) ^ ^ 

Thus ^ will be positive only if a; ^ [a;o, 0] where xq ~ —0.119 is the single real root of the qubic 

polynomial 64 x'^ — 44x2 + 36x + 5. p^j. ^ ^ ^^^^ gj always have A < ; for x ^ [xq, 0] we have 

A > 0, a decreasing function of x which varies from +oo to 5a2/64 . 

★ GMG: Combining the two theories introduce three constants that are related to the 
components of the structure tensor density (|3.4[) by 



5 (a2 + 4a6 + 4c2)^ 

A = i I A op.) 

8(a2(3a2 -40afe- 11252) -8(9a2 + 20a6 - 1662)c2 -80c4) ^' ' 

_ a^{-3a^ +A0ab + 112b^) +8{9a^ + 20ab - 16 b^)c^ + 80 c'^ 
^ ^ 16(a3 + 2a26-4a62_85c2) ^ ' ' 

a2(-3a2+40a6 + 11262) +8(9a2 + 20afe _ i6 62)c2 + 80c'* 
^~ 8(a2 +4a6 + 4c2) ^' ^ 

Inverting this system as well as discussing in general the positivity of is not very illuminating. 

For this last purpose we have plotted the region in the (x = b/a , y = c/a) plane of the sign of ^ 

on fig. [l] Let us also remark that we always have A^ = (5/64)(a2 + 4ab + Ac^)^ ^ 0. 

However, using the special values of the structure tensor density we can easily see that there is no 

solution of type O which allows to have ^ > 0; while those of type Dt and Dg are both possible. As 

for a structure tensor density of type /, we have as many parameters in the solution as the physical 

constants of the problem. Thus we may expect that for some range of values a finite number of 

solutions are always defined. Of course, when the number of geometrical parameters will be less 

than three, the solutions, if any, will exist only for special values of the physical constants. 
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Figure 1: Graphical representation of the sign of the GMG coupling constant A and as function of the rescaled 
parameters x = b/a and y = c/a in the framework of unimodular group spaces of type I. We have A = oo, § = 0, = 
on the solid curve; § = oo, A = on the dash-dotted curve; = oo on the dotted curve. On each region of the x, y plane 
delimited by these curves we have indicated the sign of A (equal to the one of ^ ) and the sign of fi. The left— hand side of 
the figure is a blow-up of the central part of the right-hand side plot. 



Type // 

★ TMG: Three different solutions occur. One, of Petrov type iV, with negative cosmological 
constant 

A = -^ , a = -^ , ; (4.29) 

9 ' 3 ' 3 ' ^ ' 

and two, of Petrov type //, with vanishing cosmological constant 

A = , a = , 6 = ; (4.30) 

A = , a = -^ , . (4.31) 

Let us mention that trivial flat solutions also occurs during the discussion of the field equations, 
for instance we meet the solution with a = 6 = 0, corresponding to the flat space |^ ; we shall not 
insist anymore on such solutions. 

★ NMG: We also obtain three types of solutions; two of Petrov type N: 

6 = , = ' ^ = (4.32) 

^This emphasises the obvious fact that a Bianchi / is a solution, always corresponds to a flat space, flat space may also 
appear as a Bianchi // model. 
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and a third one of Petrov type II : 

6=-(l±^/5)a , (61 T 1975) ^^^139^1^ 

\ J y 479 229441 ^ ^ ^ 

Let us noticed that the cosmological constant as well as the coupling constant are independent of 
u, but not the curvature of the spaces. 

★ GMG: We obtain : 

56(4a + 6)3 
8(-112a2-40a6 + 362) ' 

6(-112a2-40a6 + 362) 
^ = 16 {-4a- + 2 ab + b^) ' ^'-''^ 

6(ll2a2 + 40a6-362) 

where = (5/32)6^(40 + bf > 0. To have ^ > 0, 6 must be in the intervals ]o(20 + 
4\/46)/3, 0[U ]a (20 - 4\/46) /3, -4 a[. These solutions always are of type II. 
Solutions of type A'' also occur in the cases where 6 = 

A = , ^=^^ , (4.38) 
2a + /U 

or when b = —a 

270M + 3a . _ 17a^A^ 

^~ 272m ' ^-4{3a + 2^i) ' ^^'"^^^ 

Type III 

■k TMG: The condition /x ^ is incompatible with the field equations. 

★ NMG: There is no solution with ^ 7^ 0. 

★ GMG: There is no solution with ^ > but there is a special one of Petrov type /// with 

80 . 4000 „ , 800 „ 

a' = -7^IJ' , A = -— — u^^O , £=-;^7^M <0 . (4.40) 
69^ 12167'^ ' ^ 207^^ ^ ' 

Type IV 

★ TMG: Taking into account the condition 4i/^ > a^, we obtain as only (real) solution 

a=^Ti^ , ^=f±^ ' ^ = (4-41) 
of Petrov type Ic, and subject to the condition : u [tm/2, • 

★ NMG: Once A and ^ are expressed in terms of a, b and v, the remaining field equation 



are satisfied if 

= h(a - h) or = 

4a 
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9 ,/ ,x 9 (a^-6^)(a-6) 
i/2 = ;,(a - 6) or = ^ ^ . (4.42) 



It is easy to verify that the first is disallowed by the condition Ai/'^ > a^. Whereas for the second 
one, we may parametrise again the solutions as follows : 



b = xa , (4.43) 
iy'^^{l^x'){l + x) , (4.44) 
A.^, - , (4.45) 

^2 

C = — (8 + lla; + 18x2-5x^) . (4.46) 
8 

The condition > a? implies that x > [1 + a/5)/2 or > a; > (1 — a/5)/2 while the positivity of 
^ requires x < Xq where Xq is the single real root of the qubic polynomial 8 + 11 a; + 18 a;^ — 5 x'^, i.e. 



Xq = (^18+ v/l2987 - 60^14370 + ^3(4329 + 20/L4370)j /lb, Xq « 4.21, and insures that A > 0. 
This solution always is of Petrov type /c, both with respect to the Ricci and the Cotton- York 
tensors. 

★ GMG: The generic solution is of type Ic with 

b[{a-bY -Ay-'f 

~ "8 [(a - 6)2 (3a2 + 2Qah + 3&2) + 8z^2(a _ 9;,)(a _ _ soj/^] ' 

_ [a-bf (3a2 + 2Qab + 362) _^ g^2 _ 95) (a - 6) - 
^ " 16[(a-6)2(a + 6) -4aiy2] ' 



(a - 6)2 (3a2 + 26a6 + 362) ^ g^2(^ _ g^^^^ _ _ gQ^4 
8 [(a - 6)2 - 4iy2] 



(4.49) 



where A^ = (5/64) ((a - 6)2 - Av'^f ^ 0. We have plot on fig. [2] the zero and singular curves of 
A, /Lt and ^ on the {x — ajv, y — bjv) plane. We also find a Petrov type solution 

c 2 1 

a = 6 = , = — ^ > , with : A = - ^ . (4.50) 

2/i 5 5 



• Non-unimodular Lie Algebra 
Type T 

* TMG: We obtain two solutions, only defined for positive cosmological constant. The 
first one is of Petrov O and locally a dS-}, space, 

6 = , fc2 = A , (4.51) 
while the second one is of Petrov Type 

«=f , ^"-\^-h.^^ ' ^'^i^+l^^'' A>A^V27. (4.52) 
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Figure 2: Graphical representation of the sign of the GMG coupling constant A and as function of the rescalcd 
parameters x = ajv and y = hjv in the framework of unimodular group spaces of type lY . We have A = oo, ^ = 0, = 
on the solid curve; § = oo, A = on the two dotted straight lines; = oo on the dash-dotted curve. On each region of the 
X, y plane delimited by these curves we have indicated the sign of A (equal to the one of § ) and the sign of fi. 



■k NMG: Three different solutions are plausible. The first one is of Petrov type /r (or 
Petrov type Dg when ^ = 2 A) 



a = 



62 ^ 1(3^ +A) 



1 



fc^ = -(5A-e; 



(4.53) 



which of course requires that A > ^/5 . 
The second one is 

6 = , fc2 ^ 2 ± (e - A)) . (4.54) 

Note that for A > 0, has solution for the plus sign with ^ ^]0, A[, whereas for A < both signs 
in k'^ are accepted with ^ > 0. 

The third one, of Petrov type Dg : 

2_ - vaise +42A) 



b — a 



42 



e . vai5e+42A)-4e 



that are real only when 42A>f ^2A>0. 

* GMG: Generically we obtain solutions of Petrov type /r 

3fc4 - 1462 p ^ iQ^2 ^2 _ 



A 

e = 



2(fc2 + 8a2-5 62) 
fc2 + 8a2-562 

-fc2 -8a2 + 562 
2 



(4.56) 

(4.57) 
(4.58) 



The singular and zero curves of A, /i and ^ in the {x — a/k, y — b/k) plane are depicted on fig. |3] 
Let us mention that in the region where ^ > (the interior of the hyperbola) A > 0. 
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Figure 3: Graphical representation, in the framework of non-unimodular group spaces of type T, of the zero and singular 
curves of the cosmological constant A and of the coupling constants § and as function of the rescaled parameters x = a/k 
and y = b/k ( A = oo but 5=0 and = on the (dashed) hyperbola, A = on the black fourth order algebraic (solid) 
curve, /i = oo on the y axis ). For each region of the plane delimited by these curves, we indicate the sign of A, ^ and 5. 



There also are solutions of Petrov type D^., with — ti^ ~ a^, fi arbitrary and 



A 



2/z(17a2 +462) 



17a2 +452 



(4.59) 
(4.60) 



2{fi-3a) 

We plot on fig. |4]the zero curve of A and the singular straight line of ^ in the [x ~ a/ ^, y = 6//i) 
plane. Finally solutions of Petrov type O are obtained for arbitrary values of the coupling constant 



(+,+) 
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\(+r) 
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Figure 4: Graphical representation of the singular line {x = 1/3) of § and the zero curve of the cosmological constant A 
of the Petrov type Ds solution of GMG field equations, in the framework of non-unimodular group spaces of type T. The 
zero curve of A is a fifth order algebraic curve, admitting the vertical asymptote x = 1 and the two oblique asymptotes 
y = ±x/2. 



jjL and ^. They are dS^, (or flat) spaces but with 



A = fc2 (4 
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Type SI 

-k TMG: Two non trivial solutions occur. The first one, is Petrov O and locally an AdSs 
space with A < , 

b=-a , fc2 = -A . (4.62) 

The second one, with > A > , is of Petrov type D 

-2ii±J27A + 5ii^ , 2ii±JWX+YiJ? ,2 , /U^-27A 

a= 1 0= — , k =—ao= — 

6 ' 6 ' 36 

(4.63) 

More precisely, with the upper sign (+), the solution is of Petrov type is Dg {resp. Dt) for > 
{resp. /X < 0) ; with the lower sign (— ), it is the converse : Dt {resp. Ds) for > {resp. /x < 0). 



★ NMG: Here we obtain ^ = |[5 (a + 6)^ + 16 (a^ + 6^) + 36 fc^] > and the remaining 
field equations are satisfied in three cases. 

If 

a = -b , fc2 = 2(-C ±v/e(^ - A)) , (4.64) 

which requires that > A ^ ^ when is obtained with the minus sign in front of the square root 
or, if the plus sign is adopted : ^ > and A < or ^ < and A ^ ^ . It is of Petrov type O, 
corresponding to an AdS^ space for positive values of ^ . 

We also obtain fc^ = —ab, in which case, parametrising the solution as before by 6 = a; a, we have 

a2 4 

^ = — [8(a;2 + l) + 13(a;-l)2] > and A = — — (21 a;'' + 204a;3 - + 204 a; + 21) 

8 64 ^ 

(4.65) 

The positivity of fc^ requires that x < 0. So the maximum of A/^ is reached at x = or a: = — oo 
where the ratio tends to 1/21 . In the limit a; = 0, we obtain A = ^ /21 = As the quartic 

polynomial 21a;* + 204 a;^ — 194 a;^ + 204 a; + 21 has only two (negative) real roots : xi « —0.1 
and X2 ~ —10.7, A is negative for the value of x between these two roots. It reaches it minimum 
at a; « — 1. The coupling constant ^ increases monotonically as x decreases; when x goes to — oo, 
both A and ^ diverge. This solution is of Petrov type Dt if a; e] — 1, 0[, Petrov type O when 
x = —1, and Petrov type Dg if a; e] — oo, — 1[. 
There is also a third solution of Petrov type 7c (unless a = b = 0) 

according to sign of A we demand ^ > lOA ^Oor^ >— |A^0 

★ GMG: 
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Generically we obtain a Petrov type I 

-5(a + 6)4-8 (a2 - 30a6 + b^) P + 
~ 8 {2,0? - 26a6 + 362 _ 4^2) 

-Zc? + 26a6 - 36^ + 4fc2 



16(a - 6) 



(4.67) 
(4.68) 



^ = i (-3a2 + 26a6 - 36^ + ^k^) . (4.69) 
8 

We also recover a Petrov type O, namely an AdS^ space when 

h=-a , A=-fc2^!±^ . (4.70) 

We also obtain special solutions for = —ah, in which case we parametrise the solution as before 

by 6 = x a with a; < : 



A = -a 



2 (-21 - 55a; + 14a;2 - Mx^ + 55.t4 + 2lx^) a - 2 (21 + 204x - 194a;2 + 2Mx'^ + 21a;*) ji 



16(21(1 +a;2) - 26a;) /x 



21(1 + 0;^) -26a; 

^ 4 3(l-x)a + 2/i ■ ^ ' ' 

This solution is of Petrov type Dt if a; s] — 1, 0[, Petrov type O when x = —1, and Petrov type Ds 
if a; €] — 00, — 1[. 



Type SII 

■k TMG: We easily obtain a Petrov type iV solution 

fc = ±x/=A , a= ; (4.72) 
we also recover a Petrov type O, namely an AdS^ space when 

a = -k , A = -k^ . (4.73) 

★ NMG: Here, we obtain a Petrov type N. If 

A;2(4a + 3/c)(4a + fc) , ^ 1 /„ 2 „ , , 2\ 

^ = -WT8aW' ^ = ,{8a^ + 8ak + k^) (4.74) 

the field equations are satisfied. Accordingly, imposing ^ > 0, we obtain solutions for all values of 
a and k such that k [— 4a — 2\/2|a|, — 4a + 2\/2|a|]. All these solution have negative cosmological 
constant. Let us notice that the value of the parameter v did not play any role in the parametri- 
sation of the solution, but in the curvature. For a = —k, we reobtain a Petrov type O, the AdS^ 
geometry as a solution. For a = —k/2 we also have a conformally flat geometry, solving the fleld 
equations, but for negative value of ^ = —2a?. 
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★ GMG: 

Here, we obtain a Pctrov type N 

16aV + fc^(A: + 3A*) + 2afe(fc + 8Ai) 2 Sa^ + Safc + fc^ 

2(8^^+8^^k + ¥)Ji ' ^ ~^ 2(2a + fc + /x) ' ^ 

And a Pctrov O type for a = —k, locally an AdSa space with A = — ^ fc^ . 
Type 5 J/J 

★ TMG: We obtain a solution of Petrov type D : 

a=^ , , /x2^-9A , (4.76) 

which is Petrov type {resp. Dt) iov k = v {resp. k = —u). 
■k NMG: The following set of solutions occur: 

Firstly a solution of Petrov type Jr : 

,2 £-5A 2 3£+A 
a = , e = '^^ , v^ = , (4.77) 

which imply that A < — 3^ . 

Secondly, we also have solutions of Petrov type D : 

k = ±u , = ^ (6^+^/3^(5^ +7 A)) , i^' = ^ (4^ +v/3C(5C +7A)) (4.78) 

which requires, taking into account the condition 4 z/^ > and assuming ^ > that A ^ —5 ^ /7, 
otherwise 2 A ^ < 0. 
Finally, we also obtain : 

k = ±v , «' = ^ (ee - V3^(5^ +7A)) , !^' = ^ (4e-V3C(5C+7A)) (4.79) 

which requires that > —35^ /289 > A ^ — 5^/7 . These solutions are of Petrov type Dg {resp. 
Dt) ioTc k = u {resp. k = —v). 

★ GMG: 

Here, we obtain solutions of Petrov type J 

-5a^ - 8a2 (5z/2 _ fc^) + 16 {3k^ - lAk'^u'^ - Sz/^) 
~ 8(3a2-4A;2 + 20z/2) ' ^'^■^^> 



3a2 - 4fc2 + 20t/2 
16a 



(4.81) 



^ = ^(-3a2+4(fc2-5i/2)) . (4.82) 



We also find, solutions with k = ±v of Petrov type D 

2la^ - 42a^ /x - ISOa^ + f>lQo? iiv^ + 256a v'^ - 512/x z^' 
336/ia2 - 256/iz/' 



^ = oo...„2_orfi„„2 ' (4.83) 



-21a^M + 16AU^ (484) 
^ 4(3a-2M) ' ^ ' 
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which are of Petrov type Dg {resp. Dt) for k — v {resp. k — —v). 
Type L The solutions obtained are of Petrov type //. 

★ TMG: There is one solution , defined only for special values of ^ and A (positive): 

^ = ±3\/3A , a = -- , c = ±2\/3A . (4.85) 

c 

★ NMG: Here also only one solution is obtained, for special values of ^ and A (positive): 

^=21 A , a^-- , c = ±2%/2A . (4.86) 

c 

★ GMG: There is a class of solutions for arbitrary a 

A^-l..' , .^j^c . t^-lc- . (4.87) 

Type L' The solutions obtained are flat spaces. 

★ TMG & NMG & GMG: The only solution is the flat space, obtained with 

6 = or b=-l and A = . (4.88) 
4.2 Coordinate representation of the metric 

Having at our disposal all the homogeneous solutions of MG theories in a formal way, the purpose of 
this subsection is to perform an analytic example in which we express the invariant forms that define 
the coordinate system. The latter can be achieved algebraically by determining the GL(3,M) matrix 
which transforms the canonical expression of the structure vector and tensor density into the forms we 
use and then apply the same transformation to the canonical expression of the invariant one-forms (for 
example see the table 8.2 of ref. |13j). Another possibility consists into a direct integration of the Cartan 



equations (3.2 1 defining the invariant vector fields from which we reduce the expression of their dual 
basis. For illustrative purpose we shall perform this calculation in the framework of the non-unimodular 
SII solution, which actually corresponds to the so-called pp-wave AdS [22l[T5]] metric or null warped 
AdS metric. 

The SII structure constants correspond to a group admitting an abelian two-dimensional subgroup. 



The corresponding Cartan equations are obtained from the structure tensor density (3.23) and the space 
like structure vector ka — (0,0, fc). They read : 

[Ci,6]-0 , [iiAz] = [k + v)^^ + {u-a)^2 , [^2,^3]^-ia + iy)^i + {k-iy)^2 -(4.89) 

Thanks to the presence of the abelian subgroup, the integration of these equations is immediate (see 
ref. [13] for a discussion of this problem in the framework of the Bianchi /// group) and, after a specific 
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choice of the integration constants, leads to : 



(,i^dx , £.2^dy , £,3 = dz + r{x, y)dx + s{;x, y)dy 



where the function r{x, y) and s{x, y) are hnear 



r(x, y) — {k + v)x — {a + v)y and s(x, y) = {v — a) x + (k — v) y 



We immediately obtain the dual right invariant one-form that define the metric: 



= dx — rlx^y) dz , 0^ = dy — s{xTy) Az , 9'^ — dz 



and the metric 



{dx — r{x,y) dz) + {dy — s{x,y) dz) dz^ 



(4.90) 



(4.91) 



(4.92) 



(4.93) 



that solve the massive gravity equations when the algebraic conditions (4.721, (4.74 1, (4.75) are satisfied 



The metric (4.93) has an explicit the one-parameter (denoted A3) isometry subgroup 



z H- ^ z + A3 



but hides the two-parameter abelian isometry subgroup given by 



X ^ X + p{z) 



y^y + q{z) 



z ^ z 



(4.94) 



(4.95) 



where p{z) and q{z) are the solutions (depending on two arbitrary constants denoted hereafter Ai and 
A2 : the group parameters) of the differential system : 



whose solution reads 



p'{z) = {k + v) p{z) -{a + v) q{z) 
q' (z) = {ly - a) p{z) + {k - ly) q{z) 

p{z) — Ai e'"'^ (a cosh(az) + v sinh(az)) — A2 e*^^ (i^ + a) sinh(az) 

:= XiPi{z) + \2P2{z) 
q{z) — Ai e^^ {v — a) sinh(az) + A2 e*^^ (a cosh(oz) — v sinh(az)) 

:= Xi qi{z) + \2 q2{z) 



(4.96) 
(4.97) 

(4.98) 
(4.99) 
(4.100) 
(4.101) 



Now it is immediate to write the most general left invariant vector, a Killing vector of the metric (4.93) 



depending on the three parameters Ai,A2 and A3 introduced respectively in eqs (4.99) or (4.101) and 

i = A3 a, + p{z) a, + q{z) dy . (4.102) 



(4.94) : 
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q2 



To make explicit the action of the abehan two-dimensional subgroup we have to use the group parameters 
are coordinates. Assuming a 7^ 0, we introduce new coordinates X and Y defined by 

x=pi{z)X +p2{z)Y , (4.103) 
y^qi{z)X + q2{z)Y , (4.104) 

and we obtain 

=Pi{z)dX +p2iz)dY 

{e'''{i^ + a){dX + dY)+e-''ma~iy)dX + {a + i')dY)} , (4.105) 
= qi{z) dX + q2{z) dY 
= {e''^iy-a){dX -dY)+e-''ma-iy) dX + ia + iy) dY)} . (4.106) 

A last coordinate transformation : 

X^'L^^l±ny±^ , Y^ '^-- + k'i'^-a)u ,^_ii,(c) , (4.107) 

2k^a^/\a\ 2k^a^y\a\ k 

where a :— — sgn(a) i^, provides the usual expression of the null warped AdS metric |15) : 

1 dwdw + dC^ du^ 
d^=k2 ^^+^^2RfcTT) • (4.108) 

Obviously the coordinate system (u,v,C) (or equivalently (x,y,z)) defines a local chart but does not 
cover the all group manifold. 

4.3 Non simply transitive groups 

The Kantowsky-Saks metric, who is the metric of an homogeneous space M x S*^, on which acts (multi- 
transitively) a 4-parameter isometry group that does not contain any 3-parameter subgroup is 

ds2 ^ -dt^ + i?2 (^^02 _^ gj^2 j^2) _ {4.109) 

This geometry is conformally flat, with a traceless Ricci tensor of Petrov type Dt, whereas the Cotton- 
York tensor is Petrov type O and the Kf^^, tensor is of Petrov type Df. Thus, it can utmost satisfy the 
NMG (and so the GMG) field equations. It turns out that they are satisfied for A = , ^ ~ ■ Let 
us emphasise that the hyperbolic or the fiat version of this metric (sin is replaced by sinh 9 or hy 9 
respectively) admits a 3-parameter transitive isometry group. Thus, they have been already considered 
in the previous subsections, where indeed we find that the hyperbolic metric does not satisfy the TMG 
field equations but the equations of NMG (and GMG) are satisfied for A = ' ^ ~ ~W' 
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5 Vanishing scalar invariant geometries 



All the homogeneous geometries considered here above, share the common property, that all scalar geo- 
metrical objects which we can build are constants on the spacetime. Notably, in the frame of Lorentzian 
geometries there exist spaces (VSI spaces) which are not locally homogeneous but have vanishing scalar 
invariants built out of the curvature tensor. These geometries are a subclass of the Kundt geometries 
considered in [ID] and explicitly known in three dimension |24| . If we exclude flat space, there are two 
possibles such metric^ (labelled Al and Bl in ref. [23]) that all reads as follows 



ds^ = -2du 



dv + 2'^('"' ''^j •'') "^"^ + "^j 2;) da; 



dx^ (5.1) 



with special expressions of the function F{u,v,x) and W{u,v,x). All these metric leads to a curvature 
(Ricci) tensor that can be expressed in a null frame {l,n,m} such that l^Ua = — 1 = —m"ma as 

Rafi=(l>lJi3 or Ral3 =1pl{a'mp) (5.2) 

with 7^ and ^ 7^ 0. The first one is of Petrov type iV , the second one of Petrov type ///. Of 



course the functions appearing in the metric (5.11 can be modified by coordinate transformations which 
preserves its form: 

V 

u^U[u\ , V ^ -■ \- F[u, x] , X 1-^ X + Q[u] . (5-3) 

U[u] 

Using these transformations it is easy to check that the induced transformations on the functions {F, W) 
read 

F ^ Fu^ + 2^j-u + wgu -i'^^ -g^ , w^w + F'-g , (5.4) 

where {F, W) are expressed in terms of the (u,v,x) coordinates via their {u,v,x) dependence given by 



eqs (5.3 1. Dot and prime denote partial derivatives with respect to the new coordinates u and x. In 
what follows we briefly describe the resolution of the massive gravity field equations in case of the VSI 
geometries. Let us notice that for consistency we have to assume the cosmological constant A = 0, as 
all scalar invariants of the metric considered vanish. 

We have solved the equations of TMG, NMG and GMG for type Al and Bl. It turns out that the only 
not trivial equations are those corresponding the [u, x] and [u, u] components. Our strategy is to first 
consider the simplest, i.e. [u,x] equation and then plug its solution into the [u,u\ one and solve it. 

• Type Al Here the metric components are a priori of the form 

F{u,v,x) — V fi{u,x) + fo{u,x) , W{u,v,x) = wq{u,x) . (5-5) 



''In ref. 1241 four non flat expressions of the metric are displayed, but the last two (denoted as Dl&Fl) are special cases 
of the first ones. 
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Before proceeding to the equations of motion we shall first discuss the gauge fixing which are 



allowed by the coordinate transformations (5.3). At first we note that we can always eliminate 
wo{u, x) with an appropriate choice of T 

F{u,x) — — i wq{u,x)Ax + Q{u)x + %{u) ; (5.6) 



thus we consistently assume that wq{UtX) = 0. Using the latter, it easy to check that the residual 
coordinate transformations induce the transformations 

h=hU^2H- , (5.7) 
LA 

h = hu'' + (^xg + h) fiU'' + 2(xg + h) u-g^ 

Here after we shall present the solution of the field equations in a coordinate system fixed by 
eliminating as many as possible gauge functions in the expression of the metric. 

-k TMG The only non trivial field equations are [u,x\ and The first one gives 

dlfi{u,x) + fida;fi{u,x)=0 , (5.8) 

whose general solution reads 

h{u,x)^co{u) + c,{u)e'>'^ . (5.9) 



Using (5.7) we can eliminate cq{u) and set (locally) Ci(u) to one by choosing 2U/U = cq{u) 



Ci{u) e U = I, i.e. in Eq. (5.7) the u ~ U{u) coordinate transformation as the inverse transfo- 
rmation of 



u = 



/"e-3 /=o(«)du ^ = const. (5.10) 



and with an appropriate sign for c^, 



g{u) = - U / co{u)Au + \n[c^ci{U{u))] ) . (5.11) 



Let us notice that if ci(u) = the metric is flat. Using the latter solution, the uu equation reduces 
to 



whose general solution reads 



dlUu,x)+fidlUu,x) = -^e~^^^- . (5.12) 



fo{u,x)^qQ{u)+qi{u)x + q2{u)e ^"^ + , (5.13) 

Last but not least, we can also eliminate q2{u) with an appropriate choice of 'H{u), i.e. ^{u) ~ 
-92(")- 
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★ NMG 

Working similarly as in the TMG we obtain 



/i(M,a;) = e-^^ + C2(w)e^^ (5.14) 

and 

foiu,x) = qo{u) + q,iu)x + q^iu) e'^ ^ + ^ d^{u) ^ ^ 

6^ 6^ 

* GMG 

Working similarly as in the TMG we obtain 

;^(,,,) f e-^ + c.(.)e^-, where A, . ^W|±lZ) , ^^ ^^^ 

\ C2(w)a;e-2^^ , if ^ = -V 

and 

/o(w,x) = go(u) + ?'o('«,x) + (91 (u) + ri(u, a;)) a; + r2(u,x)e^-^ + {qziu) + r3(u, a;))e'*'+^ , 



whereas for ^ = — 

fo{u,x) = qo{u) + ro(u, x) + {qi{u) + ri{u,x))x + {r2{u,x) + (qaiu) + r3{u,x))x}e' 



-2/j, X 



I + j(u,a;) . {1 - ^ix)e^''^' j{u,x) 



e^^^ j{u, x) 



I ( \ ^ jy^, x 

r^{u,x) = — ^ , (5.18) 



where 

j{u,x) (^idji{u,x))d, + h{u,x) {dl-^d,^ ~ ^d,d^ + 2d^dlj f,{u,x) . (5.19) 

Type Bl The metric components are given by 

F{u,v,x) = +v fi{u,x) + fo{u,x) , W{v,u,x)=^ \-wo{u,x) . (5.20) 

X^ X 



Applying the coordinate transformation (5.3 1 in our case, we find that the coordinate transfo- 
rmation which preserves the form of W is the v, for which we can always eliminate wo{u,x) as 
follows 

- ^] F{u, x) - -wo(u, x) =^ F{u, x) = -x^([ dx' - niu)) , (5.21) 
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thus we can set consistently wo{u^x) = 0. Using the latter, it easy to check that the residual 
coordinate transformations induce the transformations 

F^-t,+vfi + fo , (5.22) 
fi=h-2H , /o = /o + i'H(/i-l)+2i2H 

and as in Type Al we shall present the solutions of the field equations in a coordinate system 
fixed by eliminating as many as possible gauge functions in the expression of the metric. 

-k TMG The only non trivial field equations are [u,x] and The first one gives 

,2 / I 



+ (^Ai + - j hiu, x) = , (5.23) 
whose solution reads 

fi{u,x)^co{u) + c{u)Ei{-fix) , (5.24) 

where Ei{z) = ~VV J^^e~'^/tdt, denotes the exponential integral function. Using (5.22) we 



can eliminate for example co{u) by choosing co(u) = 2'H{u). Using the latter equation, the uu 
equation reads 



'„2/ 



jiu,x) — {c^u)Ei{-nx) + 2c\u)) , (5.25) 
whose solution can be found using the method of variation of parameters 

/o(u,x) = X {qo{u) + ro{u,x) + (gi(u) + ri{u,x))x+ (g2(u) + r2iu,x)) e~^^) , (5.26) 

{1 + Hx)j{u,x) , j{u,x) , ef'^'jiu^x) 
ro{u,x) = ^ , r-^{u,x) = , r2U,x) = ^ . 

/i^X fix fl'^ X 

★ NMG 

Working similarly as in the TMG we obtain 

fi{u,x) ^ ci{u) Ei{-x^/^) + C2{u) Et{xy/^) (5.27) 

and 

fo{u,x) = X (^qo{u) + ro{u,x) + {qi{u) + ri{u,x)) x + {q2{u) + r2{u,x)) e"^^ + {q-i{u) + r3(u,a 

// N ii.u,x) , j{u,x) , e±^^j(u,x) 
r^(^u,x) = —— , r,{u,x) = — , r^^^{u,x)=T ^^2 

★ GMG 

Working similarly as in the TMG we obtain 



/i(M,a;) = I "^i*^"^-^**^^"^^ "''^2(w)-Ei(A+a;) , where X± = ^ , ^5 29) 

\ ci{u)e-^'''' + C2{u)Ei{~2fix) , if ^ ^ - 
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and 



fo{u,x) ^ X {qo{u) + raiu.x) + {qi{u) + ri{u,x))x + {q2{u) + r2{u,x)) e^^ + (qsiu) + r3(u, x)) 6-^+ ^) , 

f, s 1 + fix i{u,x) , e'^'^T j{u,x) 

ro[u,x) = ^j^j{u,x) , r^{u,x)^ , r^isiu^x) = ^ {X+ - X )x ' ^ 



whereas for ^ = — 

fo{u,x) = X (qaiu) + rQ{u,x) + + riiu,x))x + {q2{u) + r2{u,x) + (qsiu) + r3(M, a;))a;} e~^^^) , 

1 + fix , j{u,x) , {1 - fixje^^f" j{u,x) 

ro[u,x) = — J[u,x), r-^{u,x) = —^, r2[u,x)- 



e^f^^ j{u, x) 
4/i2 X 



r^{u,x) = —-^^ , (5.31) 



where 



jiu,x) ( id,h{u,x))d., + hiu,x) ( dl-^9,j ~ ^dA + 2dud'', ] h{u,x) . (5.32) 



6 Conclusion 



To summarise, we have obtain all homogeneous spaces, solutions of the TMG, NMG and GMG field 
equations (with cosmological constant), at least formally. We classified them according to canonical 
representations of the structure constants obtained by Lorentz transformations. To obtain the explicit 
expression of the metric, we still have to solve an elementary algebraic problem, as this was explicitly 
demonstrated in subsection |4.2| for the type SII, which consists to find the linear transformation that 
maps the expressions of the structure tensor density and the structure vector obtained on their usual 
canonical expression, i.e. express the frame in the coordinate basis {9°'[x^\) as linear combination of the 
canonical one or to directly integrate the expressions of the right invariant forms. 

We also have determined the Petrov types of the traceless Ricci tensor, the Cotton- York tensor and 
the traceless K^^, tensor of all Lorentzian 3-dimensional homogeneous geometries, which proves to be a 
useful tool to recognize equivalent solutions: for instance solutions of TMG, which are of Petrov type D, 
are biaxially squashed AdS^ geometries |9j. 

In brief, we found that those of type / (which corresponds to all unimodular Bianchi types) allow all 
values of the cosmological constant; for types // (which could lead to the Bianchi types //, VIq, VIIq 
and VIII; for types /// and IV (which correspond to Bianchi type VIq and VIII), in type /// only the 
GMG has a solution of Petrov type /// with negative cosmological constant. In case of homogeneous 
spaces corresponding to non-unimodular algebras, the Lorentzian classification of the algebra leads to 
solutions of types T and SI (which correspond to all non-unimodular types Bianchi types ///, IV, 
Vlh and Vllh) which allow both signs of the cosmological constant; for type 5// (which correspond 
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to Bianchi types IV, III and VIk) with only negative values of the cosmological constant; for type 
SHI (which correspond to Bianchi types /// and Vlh) which allow both signs of the cosmological 
constant; for type L (which correspond to all non-unimodular Bianchi types, depending on the values 
of ac) with positive cosmological constant for TMG & NMG and negative for the GMG; for type L' 
(which correspond to Bianchi type V for h — Q and to ///, Vlh for b 7^ 0.) 

In addition, we have also obtain the solutions of TMG, NMG and GMG field equation for VSI 
geometries (which imply a vanishing cosmological constant). The results are found to be of Petrov type 
///. They all contain, after having fixed the coordinate system, several arbitrary functions of a lightlike 
coordinate, reflecting the physical degrees of freedom of these solutions. 

The above results provide a classification of all the solutions of massive gravity theories on CSI 
geometries. Of course not all of them can be considered as classical background geometries. For instance, 
it is well known that on Bianchi IX solutions, due to the compactness of the space, no global causal 
structure could be defined. It would be very interesting and instructive to study the perturbative stability 
of these backgrounds, and identify all the physically relevant configurations. It would be newsworthy 
to check if the Petrov type D solutions of NMG and GMG are biaxially squashed ^4^5*3, like as in the 
TMG [1]. 
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A Petrov classification of homogeneous spaces 

In this appendix we provide the Petrov classification oi Sap , Ca p , Ka p tensors, for homogeneous spaces 
of types A and B, corresponding to unimodular and non-unimodular Lie algebras. 

• Tables 1 and 2 correspond to traceless Ricci tensor for types A and B respectively. 

• Tables 3 and 4 correspond to Cotton- York tensor for types A and B respectively. 

• Tables 5 and 6 correspond to Ka/i tensor for types A and B respectively. 
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Finally, we give the expressions of the discriminants of p tensor for all the cases 

= -TT7T^((a + bf - c^fc^ (Sa^ - 2a% + AOab'^ + 646^ + 4(5a + 2mc^f x 
" 110592 

X (21a^ + 4ac(-66 + 5c) + 0^(46 + 26c) + 8c (86^ - 86c + bc^)f x 

X (21a^ + a^(46- 26c)+4ac(66 + 5c)-8c(86^ + 86c + 5c^))^ < , (A.l) 

= , (A.2) 

\f = , (A.3) 

A'/ = —4- Ui^^ - a") (-ab + b^ + v^f {-21a!' + lha% + ab'' + 56^ + 4(13a - bb)y^f x 
442368 

(I05o^ + IQaH + 381a^62 - leOa^fe^ - 450^6"^ + 84a6^ - 4416^ - 4 (iTla" + 560^6 + 4860^62 - imab^- 
-416^) + 16 (91a^ - 2Qab + 1056^) z/^ - imOv^f > , (A.4) 
A| = -^6^ _ ^2 ^ ^2^2 |^^g^3 ^ ^g^^2^2 ^ ^^^^^4 _ 740^6^ + 256^) + (49a2 _ m^) + k^) x 

X (64a^ + (-56^ + ^2) V 160^(62 + 3fc2))^ <0 , (A.5) 

= 5^ — (a + ft)** (a6 + k'^f (105a'' - 1560^6 + 5020^6^ - 15606^ + 1056'' - 8 (iQa^ - 58a6 + 196^) 

^ 442368 ^ ' ^ ' ^ ' 

+16fc^)^ [{a - bf (21a^ + 10a6 + 216^)^ - 4 (63a'' - 348a^6 + 9700^6^ - 348a6^ + 636^) 

16(39a2-118a6 + 3962)fc''-64fc'^) , (A.6) 

Af ^ = , (A.7) 
Af " = {Av^ - a'f (fc2 - v^f (l05a4 + 16 {k^ - ^v^f - 8a^ {l9k^ + SSi/^))" x 

X (AAla^ - 84a* {Sk'^ + 26i/^) - 64 {k^ - bkv'^f + IQa^ (39fc* - 24A;^z/^ + 169;/*)) , (A.8) 

A| = . (A.9) 
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